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Abstract 

This paper investigates a background charge one Skyrme field chirally 
O I coupled to light fermions on the 3-sphere. The Dirac equation for the sys- 

tem commutes with a generalised angular momentum or grand spin. It can 
be solved explicitly for a Skyrme configuration given by the hedgehog form. 
The energy spectrum and degeneracies are derived for all values of the grand 
Tj- I spin. Solutions for non-zero grand spin are each characterised by a set of four 

polynomials. The paper also discusses the energy of the Dirac sea using zeta 
function regularization. 



> 



od 

a^ : 1 Introduction 

o 

The Skyrme model is a nonlinear SU (2) field theory which gives a good description 
^ . of atomic nuclei and their low energy interactions [1]. In addition to the funda- 

mental pion excitations the theory also has topological soliton solutions, known 
as Skyrmions. These are labelled by a topological charge or generalised winding 
number B, which can be interpreted as the baryon number of the configuration. 
On quantization, Skyrmions are found to describe nuclei, A-resonance [2] and also 
bound states of nuclei, see [SI HJ El El [TJ [H] for the quantization of multi-Skyrmions 
and [HI [TOl [11] for recent quantitative predictions of the Skyrme model. It is well 
known that Skyrmions can be quantized as fermions [121 [13]. Therefore, when the 
Skyrme field is coupled to a fermion field, there are two different ways of describing 
fermions in the same model. The fermion field can then be thought of as light quarks 
in the presence of atomic nuclei, [Hj. In the presence of a Skyrme field the energy 
spectrum of the Dirac operator shows a curious behaviour, namely, a mode crosses 
from the positive to the negative spectrum as the coupling constant is increased, 
[15j . In a very similar model, Kahana and Ripka calculate the baryon density in the 
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one-loop approximation [16] and the energy of the Dirac sea quarks [T7]. Recently, 
these calculations have been extended to multi-Skyrmions, [T5| 

Static field configurations in the original Skyrme model in fiat space are given 
by maps — SU(2). By using the boundary condition for the Skyrme field to 
unify the domain of such a map with infinity we make the domain compact and 
equivalent to the 3-sphere S^. If we also consider that is the group manifold 
of SU{2), we can see that the field configurations are topologically equivalent to 
maps and, because of this, the model can be generalized to the base space 

being a sphere of radius L, [20l [21]. In the limit L oo, the original model is 
recovered. The advantage of working on is that the Bogomoln'yi equation can 
be solved for B = 1, and the solution is given by the identity map |20]. This enables 
us to calculate the energy spectrum and the corresponding fermion wave functions 
explicitly. 

In [22] a system of light fermions, on S^, coupled to a spherically symmetric 
background Skyrme field was studied for grand spin G = 0. In this paper we consider 
the general case where the grand spin also takes positive integer values. The Dirac 
equation on is derived in section [2] through the use of stereographic coordinates. 
In section flTH the solution of the Dirac equation for G = is reviewed. In sectional 
the correct ansatz for the spin-isospin spinor for general G is deduced using parity 
arguments. We then present the general solution. Plots of energy against fermion- 
Skyrmion coupling constant are also given. In section [3^ we discuss the degeneracy 
of energy eigenvalues. In section HJ we address the problem of calculating the energy 
of the Dirac sea using zeta function regularization. We end with a conclusion. 



2 The Dirac equation on 



Following [22|, we now recall the Dirac equation on a 3-sphere of radius L = 1. 
Consider the stereographic projection from the north pole to the plane through 
the equator. Let be embedded in with coordinates {xi,X2,X3,w). As a result 
of projecting from onto the equatorial plane, points of can be labelled with 
coordinates Xi. The chart is defined everywhere apart from the projection point, 
A^. The coordinates Xi can be written in terms of coordinates as 



1 — w 

We define = Xl + X| + Xj. Then the metric can be written as 

/ 4 4 4 \ 

^?Mx53 = diag ^1, -^YTW " (T + W " (T + W J ■ 

We now choose the non-coordinate basis 

ea = e^^dx,. (3) 

It is convenient to choose diagonal vierbeins Cq,'^, such that 

eo" = 1, Ci' = -— , (4) 
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where all other components vanish. With our choice of vierbeins, we can calcu- 
late the matrix valued connection 1-form uJap- The 1-form uJap satisfies the metric 
compatibility condition uja/s = —^^a, and the torsion- free condition 

dO" + u^p A = 0, (5) 

where 6'" = d^^dX^ is the dual basis of Cq. After a short calculation we find 

-{X'^dXl^ - X^dX'') a,/? = 1,2,3, 



otherwise. 

The spin connection Vt^ can now be expressed as 

n^dX^' = (7) 

where = | [7a, 7/3] and the components of the commutator are the standard 
gamma-matrices, satisfying {70,, 7/3} = Irjap. We work with the following represen- 
tation of gamma-matrices 

^ V ; ' ^ ^ -a, J ' ^ ^ I2 

because we will be working with parity eigenfunctions. Here cXj denotes the set of 
three Pauli matrices, defined by 

For massless fermions in curved space-time, the Lagrangian is 

^fermion = ^ («7°ea'' (<9« + (10) 

With our choice of coordinates and vierbeins, we obtain 

Cferrmon = t) (^t^'dt - (^^^"^x, - X^ ^ t). (11) 

In this paper, we investigate fermions coupled to Skyrmions on S^. We consider a 
background B = 1 Skyrme field coupled to the fermions. The full Lagrangian C 
is the sum of the fermion Lagrangian Cfermion, the Skyrmion Lagrangian Cskyrmion 
and the interaction Lagrangian Cint- We consider fermions in the background of a 
static Skyrme field and neglect the backreaction. Therefore, we no longer discuss the 
Skyrmion Lagrangian, and the interested reader is referred to [H]. Cm is derived 
in [23], namely 

Cint = -gi'io- + i-f^T ■ 7z)lp, (12) 

where U = cr + ir • tt is a parametrization of the Skyrme field and g is the coupling 
constant. ip is a. spin-isospin spinor. It is convenient to split the spinor into two 
2x2 spin-isospin matrices ipi and such that 
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Since any complex 2x2 matrix can be expressed as a linear combination of the Pauli 
matrices and the identity, it is convenient to choose these four as a basis of SU{2). 
The spin-isospin matrices can then be written as ipi = ao^l2 + iaf'ak- With this 
notation spin operators act on by left multiplication, ak'ipi, whereas the isospin 
matrices act on by right multiplication, 

= -i^i(T2(rk(^2- (14) 

In this paper, we only consider spherically symmetric Skyrmions. The B = 1 
Skyrmion on is spherically symmetric [2D], but for B > 1 this is no longer true. 
Spherically symmetric Skyrme fields are best expressed in terms of polar coordinates, 

U = exp(i/(/i)e^ ■ r), (15) 

where f{fi) is the "radial" shape function and is the unit vector in the /x direction, 
see equation (fTOl) . Using (fTTjl and (fT2|) we can write down the Dirac equations for 
fermions coupled to a spherically symmetric background Skyrmion. We obtain 

(^1% - (^-^dx^ - - 9U'') t) = 0, (16) 

where 

t/Ts = cos/(/i) + ^756^ ■ rsin/(/i). (17) 



2.1 Solutions of the Dirac Equation for G = 

The above Dirac equation (fT6l) and the ansatz %l){Xi,t) = e^^^ip^Xi) lead us to the 
time-independent Dirac equation 

^ / gcosfifi) (7 ■p + igSi.-T sin f{fi) 

\ crp- ige^ ■ r sin /(/i) -g cos /(/x) 

where 

(sin 9 cos 
sin 9 sin 
cos^ 

and 

(7 ■p = -t(e.- (7 (d. + -SiRJ^] - -^eg ■ adg - . \ . e^ ■ crOs) . (20) 

■t^ V \ ' 1— cos/^ / sin^ " " sin /i sin y V VI \ ' 

The elements of the matrix in fllSp commute with the total angular momentum 
operator G = L + S + I where L is the orbital angular momentum, S = |cr is the 
spin operator and I = is the isospin operator. 

The above equation ( |T8l) is invariant under parity P where 

Pij{X,) = 7oV^(-X.), PX,P"i = -X,. (21) 




cos ^ cos \ 




1 — sin \ 


COS 9 sin 




cos 


— sin 9 1 




^ / 



(19) 
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The G = case is treated in [22]. There the ansatz for gives rise to a system of 
two first order ODEs, which can be expressed as a second order ODE. This equation 
can be solved analytically for /(/i) = and /(/x) = /x. In [22] the following energy 
spectrum was derived for /(/i) = 0, 



E = ±^g^+{N+I) for iV = 0,l,2,... (22) 

Setting u = cos/i, the eigenf unctions Gn{u) were found to be given by Jacobi 
polynomials. The shape function /(//) = /x was also considered in [22] • This leads 
to a second order Fuchsian equation with four regular singular points, two at n = ±1, 
one at infinity and one depending on E and g. The equation could still be solved in 
terms of polynomials. The following energy spectrum was derived, 

Eo = l-g, Et = \±^n'' + 2n + {g- If for n = l,2,... (23) 

with eigenf unctions 

G'„(n) = 5];a,(n + 1)^ (24) 

j=0 



where 



i!(2j+l)!! IIV^ - 1 y ' 4 

1=1 

for i = 1,2,... and ao = 1. Here (2j + 1)!! = 1-3 (2j + 1) is the product of 

odd integers. Using E=l± ^{n + lf -2g + g'^ in the product in (125|) . aj can be 
written as 

°. = '''''''l;;f]^"r'^' n((''-')("+'+2)). (26) 



Expanding the product in 0261) we obtain 



I j / ('«+l)K2i+l)!! n{n+2) ' ' > 



3 Solutions of the Dirac equation for general G 

In the following we derive the solution of the Dirac equation (fTS!) for general G. 
As a starting point, we construct the total angular momentum operator eigenstates 
\jm)^ and \ j'm)2 in terms of angular momentum and spin states. They are expressed 
as 

= \\-\)s + IH>5' (28) 

where Yj^m is a spherical harmonic and || ± is a spin state. For general G we 
consider four eigenstates, each of which can be written in terms of and \jm)2- 
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These are 



\GM)^, = \j=G-lm = M + \)^^^ \\-\)^ 

'^|j = G-l,m = M-i>^Jli>^, (30) 



2G 



and 



|G^)m = V^W^ |j = G + I, m = M + \)^^^ I \ - 1)^ 

^^\j = G + \.m = M-\)^^\\\)^^ (31) 



2G+2 



where || ± is an isospin state. To carry out our analysis, an ansatz in terms of 
the G-eigenstates must be found. For each row of our ansatz for the spin-isospin 
spinor the states must be of the same parity. Under parity 



so that 



It follows that 



Yi,m ^ {-l)%m. (32) 

|jm),_2 - (-1)^'^^ (33) 



\GM)^ (-1)^-1 |GM)„, \GM), {-ir |GM)„ 
\GM)^ {-ir |GM)^, \GM), (-1)^+1 \GM)^ 



(34) 



We see that \GM)^^ and \GM)^ both have parity (-1)'^ and that \GM)^ and \GM)^ 
both have parity —(—1)*^. Hence the spinor 

/ h{u)\GM), + c{u)\GM)^ \ 
^ \d{u)\GM), + a{u)\GM)J^ ^'"^ 

will have overall parity (—1)*^. Clearly exchanging the upper and lower rows will 
change the parity by a factor —1. A short calculation shows that this is equivalent 
to making the transformation 

9^-9 (36) 
in the resulting equations. For general G and parity (—1)*^ we make the ansatz 



VT^VT^^G2{u) \GM)^ + ViTTlVT^^ ^Gsiu) \GM)^ 
tVTTTLVT^^G^iu) \GM)a + tVT^VT^^~^Gi{u) \GM)^ 



(37) 

where Gi{u), G2{u), G^i^u) and G4^{u) are functions of u to be found, and the 
normalization factors are chosen for later convenience. Substituting this state into 
( IT8|) . we obtain a system of four coupled first order differential equations in Gi{u), 
G2{u), Gs{u) and G^i^u). We will solve for the case /(/i) = fi. 

We need to know how the operator ■ cr acts on the G-eigenstates in (l30ll and 
( I3T1) . The results we require are 

e,.cT\GM)^ = -\GM)^, e, ■ a \GM)^ = - \GM)^, (38) 
e,-cT\GM), = -\GM),, e, ■ a \GM) , = - \GM),. (39) 
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These results can be derived by first deducing that 



e 



cr\jm)-^ = -\jm)^, (40) 
e^,-CT\jm)^ = -Ijm)^. (41) 

To obtain the relations (l40l) and (14T!) the operator ■ <t is expressed in terms of 
spherical harmonics. Formulae for products of spherical harmonics are then needed. 
The required results can be found in 0381) and 0391) then follow. 

We also need to know how the operator {—eg ■ crdg — ^j^e^ ■ crd^) acts on the 
G-eigenstates. The results are 

i-ee ■ ado - -^^e^ ■ crd^) \GM)^ = -(G - 1) |GM)^, (42) 

i-ee ■ cTde - -^^e^ ■ ad^) \GM), = -G |GM)„ (43) 

i-ee ■ CTde - -^e^ ■ crd^) \GM)^ = (G + 1) |GM)„, (44) 

i-^o ■ CTde - -^e^ ■ cTd^) \GM), = {G + 2) \GM),. (45) 

In order to derive equations ( H2l) -( H5l) we note that 

2L ■ S = S+L_ + + 2S3L3, (46) 

where 5*+ and are defined as S+ = 81+182 and 8^ = 81 — 182 and (5*1, 5*2, 5*3) are 
a set of generators of the Lie algebra of 8U (2) and are related to the Pauli matrices 
via 8i = L+, L_ and L3 are the orbital angular momentum operators. We note 
the following result 

i-ee ■ CTde - -^e^ ■ ctO^) = (e^ ■ ct){2L ■ S), (47) 
which can easily be proved by multiplying out. Then 

2L-S|jm), = (j-i)|jm),, (48) 

2L-S\jm), = -U + I)\jm),, (49) 

can be proved by considering how L+, L_ and L3 act on the spherical harmonics. 
The necessary formulae can be found in [23]. These two equations also follow from 
2L ■ S = J2 - - S2. 

It can then be seen that 

{-ee ■ CTde - -^e^ ■ crd^) \jm)^ = -(j - i) \jm)^ , (50) 

{-ee ■ CTde - ■ o-d^) Ijm)^ = (j + |) \jm)^ . (51) 

This leads to (US)- (US}. 

The operator ■ r acts on the G-eigenstates to give 

e, . T \GM)^ = -^^^^1^ \GM), - |GM), , (52) 

e, . T \GM), = -^^^^1^ \GM)^ + |GM), , (53) 

■ T |GM), = -^J^ |GM)„ - |GM), , (54) 

■ T \GM), = |GM), - ^^^^^^1^^ |GM), . (55) 



7 



These equations can be proved by expanding ■ r then multiplying this, from the 
right, by each G-eigenstate. The right hand side of each equation is then computed 
and matrix components are compared. 

Substituting ([37]) into ([IB]) and using the identities ([M]), ([SH]), ([12])-([1S]) and (13^- 
( 155]) we are led to two equations, one in \GM)^ and \GM) ^ and one in \GM)^^ and 
\GM) ^. We equate coefficients of the G-eigenstates to obtain four ODEs in 
G2{u), Gslu) and Gi{u). These ODEs are shown below. 



dG, f 9{l-u) \ 2g^G{G + l)il-u') 



(1-^.)^ - [G+-,-^^^\G,-{E + gu)G,+ ^ ^ G„ (57) 



dG, / 3 gil + u) \^ 2gVG(G + l) ^ 



These equation give the solutions for states with parity (—1)*^. Due to the symmetry 
f l36p states with parity (—1)'^"'"^ are obtained by replacing g by —g in the equations 
above. 



3.1 The energy spectrum 

In this section, we derive the energy spectrum of the time-independent Dirac equa- 
tion (ITS]) and the corresponding eigenf unctions. We discuss a useful symmetry of 
our system of equations (ISB]) - (ISI?]) and also comment on associated second order 
and fourth order equations. In order to derive the spectrum, we use the theory of 
Fuchsian differential equations, and in particular, regular singular points and their 
exponents, see [25]. Finally, we present the explicit solution. 
Under the transformation 

G2(m), Gsiu), G,{u)) ^ i-Gsi-u), -G,{-u), -G,i-u), -G^i-u)) (60) 

followed by M I— > —u, f l57p is mapped into fl59p and fl56p into fl58p . and vice versa. 
Hence the system of equations (1S^ - (|3U]) remains invariant. Eliminating Gi{u) and 
Gslu) from the system (|56]) - (|59]) results in a system of two second order equations, 
which again map into each other via the symmetry ([60]) . These second order equa- 
tions prove to be useful for deriving an ansatz for G2{u) and G^^u). Finally, we can 
also derive two fourth order ODEs, by eliminating G^^u) and G2{u), respectively, 
again related via (jBD]) . Hence, if we have solutions for G2{u) we can find solutions 
for G^i^u). We can then use our system of first order ODEs to find solutions for 
Gi{u) and G3{u). 

In order to solve our equations we first have to derive the energy eigenvalues. 
Our system of equations ([5B]) - ([S9]) has regular singular points at m = ±1 and an 
irregular singular point at infinity. For both regular singular points the exponents 
are 

(0,0,-G-l,-G-|). (61) 
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We require our solutions to be regular over the whole 3-sphere and, in particular, at 
the north and south poles, u = ±1. The solutions corresponding to the exponents 
—G — I and —G — | contain poles, so we can exclude these solutions. The fact 
that there are two exponents taking values of zero means that corresponding to 
each regular singular point is a solution with logarithmic terms and the solution is 
therefore singular. As a result, we can also exclude these solutions. The regular 
solution can therefore be expanded as a power series in 1 + m around the south pole, 
and also a power series in 1 — m around the north pole. These two expansions only 
agree for certain values of the energy E. It turns out that these energy eigenvalues 
can be calculated from the exponents at infinity of the fourth order ODE in G2{u), 
mentioned above, which arises by eliminating Gi{u), Gs{u), and 6*4 (u) from our 
system of equations flSBl) - flS^ . This equation is of Fuchsian type and has three 
regular singular points at n = ±1 and infinity. 

The solutions of Fuchsian differential equations can only have singularities at 
their singular points. According to our discussion above, we are interested in the 
solutions of (ITSl) which are non-singular. Therefore, G2{u) has to be regular at 
u = ±1, and hence on the entire complex plane. So, G2{u) is an analytic function, 
in fact an integral function, on the complex plane. As it is the solution of a Fuchsian 
differential equation, it can only have poles at infinity, and it follows that G2{u) is 
a polynomial. 

The exponents corresponding to m = oo can be found by setting u = - then 
considering z ^ 0. We obtain the exponents 

= l + G±|Vl + 4E2 + 4E- V, (62) 
Pa = 1 + G ± i + 4^2 - + 8(7 - V. (63) 

As argued above 6*2 (^) is a polynomial. Let its degree be denoted by k. Then 
the exponents at n = oo can be equated with —k. From the exponent Ps we obtain 
the following energy eigenvalues 

= -i ± v/(fc + G+TF+7 forG = l,2,..., k = 0,l,... (64) 

Efy^ is a novel feature which arises for G > only. Note that this energy is invariant 
under g —g. From pa in (!63l) we obtain another family of energy eigenvalues, 
namely, 

E^sym = -2 ± V{k + G+iy-2g + g^ for G = 0, 1, . . . , k = 0, 1, . . . , (65) 

where G and k are not both zero. This energy spectrum has already been obtained 
in [22] for the case G = 0. A slight subtlety occurs for k = and G = 0. In this 
case, only 

Eo = l-g (66) 

leads to a regular solution. The energy level (|66l) is rather special as it crosses from 
the positive spectrum to the negative spectrum as the coupling constant g is varied, 
also see [22] for further details. 

Now that we have derived the energy spectrum, we can solve the system fl56l) - 
fl59|) by first considering the fourth order ODE in G2{u). We make the ansatz that 
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G2{u) is a polynomial in 1 + m and insert this into our ODE to find the polynomial 
coefficients. The symmetry (MJh and the system of second order equations for 6*2 (^) 
and G^i^u) leads to a related expression for G4^{u). Then the solution corresponding 



to E^gy^ is found to be 



fc k 

G2{u) = Y,aj{l + uy and G4H = (-1)'= ^ aj(l - n)^'. (67) 

j=0 j=0 



The general expression for aj is 

_ /_.y f k \ (2G+fc+j+l)!(2G+l)!! {E+g~^-G){E~g+^+G) 
V V j / (2G+fc+l)!(2G+2j+l)!! 2Gik+g)+k{k+2) ' 

If we set G = 0, and hence k = n, the above formula leads us to (!27|) which is 
equivalent to the result from [22] . 
For E±„ we find 

fc fc 

G'2(u) = ^aj(l + u)^' and G'4(u) = (-1)'=+^ ^ 0^(1 - u)^'. (69) 
i=o i=o 

The general expression for aj is now 

n = r-lV+l ( ^ \ (2G+fc+j+l)!(2G+l)!! {E+g+^+GXE-g-^-G) . . 

V I j / (2G'+fc+l)!(2G+2j+l)!! 2{G+l)(g-fc)-fc2 • y'^f 

We then use equations fISUl) and (IFTIl to obtain Gi(m) and G3(m), respectively, and 
it is easy to see that Gi{u) and G^i^u) are polynomials of order k + 1. 

We can carry out a consistency check on our solutions by setting g = in 
equations fl56|) - fl59|) and manipulating the equations to obtain two uncoupled second 
order ODEs. These are both Jacobi equations and have polynomial solutions which 
can be expressed in terms of hypergeometric functions (see [2S])- For g = 0, our 
solutions are the same. 



3.2 Degeneracy of the energy spectrum 

In order to discuss the degeneracy of the energy spectrum it is convenient to intro- 
duce n = k + G, where the integer n is analogous to the principal quantum number 
arising in the quantum mechanics of the hydrogen atom. Then the energy spectrum 
for states of parity (—1)*^ is given by the two families 

Efyjn) = -l±^(n + l)^ + g^ for n = 1, 2, . . . , (71) 
E^^y^in) = l±^(n + iy-2g + g^ forn = l,2..., (72) 

and the special energy level (l66l) . 

Eo = l-g. 
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Figure [T] shows the energy spectra for different values of n. There are two different 
ways of reading figure [TJ The obvious interpretation is the energy spectrum of 
states of parity (— 1)*^ as a function of the couphng constant e M. For the second 
interpretation and in the following, we restrict our attention to g > 0. Then, the 
negative values of g corresponds to states with parity (—1)*^"^^ due to symmetry 
(!36|) . while positive values of g again correspond to states of parity (—1)'^. The 
latter interpretation is very useful for discussing the degeneracy of the spectrum. 




Figure 1: The energy E a.s a. function of the coupling constant g for Eq (solid red 
curve), E^ (n) (dotted blue curves) and Ef {n) (dashed yellow curves). 



The energy level ( !66l) only exists for G = 0. It gives rise to a positive parity state 
with energy \ —g and a negative parity state with energy \ + g- Since the degeneracy 
of a state with grand spin G is 2G + 1 these two states are non-degenerate for (7 > 0, 
and "parity doubling" occurs for (7 = 0, |14j . 

For E^gy^{n), positive and negative parity states will in general have different 
energy eigenvalues for a given value of the coupling constant g. Recall that n = k+G, 
hence G can vary from to n. Therefore, the degeneracy is 

n 

D{E^symin)) = Y,{2G + l) = {n + if. (73) 

G=0 



11 



For Efy^{n), positive and negative parity states have the same energy for a given 
value of the couphng constant g. These states only exist for G > 0, hence G varies 
from 1 to n. Therefore, the degeneracy is 

n 

D{E%^{n)) = 2 J](2G + 1) = 2n(n + 2), (74) 

G=l 

and the extra factor of 2 is due to parity. 

We now consider the case of zero coupling [g = 0) which is equivalent to massless 
fermions on S^. In this case there will clearly always be invariance under g — * —g, 
so 

Dg=o {Eo) = 2, Dg=o {Etsymin)) =2{n + l)\ Dg=o {E%^{n)) = 2n (n + 2) . (75) 



At = I the energy level Eq and £"+^(1) are degenerate, hence the degeneracy is 

D^=,{E = I) = 8. (76) 

The energy eigenvalue £" = — | only occurs for E~^y^{\), hence the degeneracy is 
again 

3 



Dg=,{E = -I) = 8. (77) 

ivmuc — iv -r ^ 

E+{N + l), hence 



The energy eigenvalue E = N + N = 1, 2, is attained by E^^y^{N), and 



Dg=o {E = N+l)=2{N + 1){N + 3) + 2{N + if = A{N + 1){N + 2). (78) 

Similarly, E = —N — | is attained by E^gy^{N + 1) and E~y^{N), hence 

Dg=o {E = -N - ^) = 2{N + 2)2 + 2A^(iV + 2) = 4{N + 1){N + 2). (79) 

After considering the factor of 2 due to isospin and another factor of 2 due to parity, 
equations flTSl) and fl79|) are consistent with the results in [27]. 

So far, we have only considered generic degeneracies and the case g = 0. This 
energy spectrum is rather special in that we can also calculate all the accidental 
degeneracies for g > 0. These degeneracies all occur for rational values of g. For 
example, the negative parity state with energy ^ + g is only degenerate with the 
states with Efy^{n) for 

g = \{n-l){n + d,) (80) 
and with the (—1)*^+^ parity states with E~gy^^{n) (changing g to —g) for 

g = \n{n + 2). (81) 

The positive parity state with energy | — is always non-degenerate for g > Q. 

( 

asym \ 



Similarly, (—1)*^ parity states of energy E^ (n) are degenerate with (— 1)*^"^^ parity 



states of energy E^^y^{h) for 

g = ln{n + 2)-ln{n + 2), (82) 
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which is positive for n > n. Finally, states with energy E^^y^{n) and states with 
energy E'^y^in) are degenerate for 

^ ^ 4(n+l)^-(l + (n + l)^-(n + l)Y 

^ 4(l + (n + l)2-(n + l)2) ' ^ > 

and a similar equation holds for E~ {n) and E~ {n). 



4 The Dirac Sea 

In this section, we briefly comment on the zeta function regularization [28]. In order 
to calculate the energy of the Dirac sea, 

CXD 

^Dirac = 5^/^(iV)^(iV), (84) 
Af=0 

where E{N) is the A^th negative energy eigenvalue, see e.g. [29], and D{N) is its 
degeneracy, we define the zeta function 

oo 

as) = Y.D{N)EiN)-^ (85) 

N=0 

The expression (|84l) is clearly divergent. However, the expression ( !85l) is convergent 
for large enough s. The Dirac sea energy (18^ is then defined by the analytic 
continuation of (l85i) to s = — 1. For example for = we have 

oo 

E,=, = -4j2{N + l){N + 2){N+lr^ (86) 

s= — l 



N=0 

Hence, the relevant zeta function is 



C,=o{s) = -4 J]((iV + -1){N+ |)-^ (87) 



N=0 

which can be rewritten as 



where Ch{s, a) is the Hurwitz zeta function defined as 

oo 

a(s,a) = J](n + a)-^ (89) 

n=0 

Evaluating (g=o{s) at s = — 1 we obtain 
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For massive fermions on corresponding to the case /(/i) = in Section \2A\ 
the energy is given by fl^2]) and the degeneracy is 

D{N) =4{N + 1){N + 2). (91) 

Hence the Dirac sea energy is given by 

oo 

Efi,)=o = -4 ^ (iV + 1) (iV + 2) ( (iV + I) % "1 (92) 

N=0 ^ 2 

Zeta functions of generahzed Epstein-Hurwitz type are of the form 

oo 

F {s; a, b^) =J2{in + af + b^) . (93) 

n=0 

Asymptotic expansions for are discussed in [301 US\- Here, we are concerned 
with a generahzation of fl93|) . namely, 



(s; a, 6') = ^(n + a)" {{n + + b^) \ (94) 

n=0 

where we assume that a > and 6 > 0. We follow [28] to derive a formula for 
F^"^\s = — |; a, 6^). We first perform a binomial expansion which is valid for b < a 
and rewrite (!94|) as a contour integral 

n=0 fc=0 I + J l-^J 

= V— /' r(s + ^)fe^^(n + a)-^^-^^+"^ TT 

n-O fj 

Recall that 

TT (-1)^ 

— -— = ^ '- + Oiz - k) for k e Z. (97) 

sm(7r^) z — k 

The contour C encloses all the non-negative poles of 1/ sin(7r2;) with anti-clockwise 
orientation and can be split into a part 

—ZQ—ioo 



-ZQ+ioo 



where < 2:0 < and a semi-circle at infinity. The latter does not contribute to the 
integral. Now, we can move the sum over n under the integral and use the definition 
of the Hurwitz zeta function 0891) to obtain 



—ZQ—ico 

^"'""'^ 2z J T{z+l)T{s)sm{7rz) 

14 



This can be evaluated by closing the contour again, and using Cauchy's theorem. 
This time the contribution of the integral over the semi-circle at infinity is non-zero. 
However, it was shown in [30] that the contribution is very small, so we neglect it 
in the following. 

From now on, we focus on the physically relevant value s = — |. The integral (l98ll 
has poles at 2; G Z due to 1/ sin(7rz). Only the non- negative poles contribute, because 
of the contour. The gamma function T{z — ^) has poles at 2; — | = 0, —1, —2, . . . 
Only the pole at z = | lies inside the contour. Finally, there is a contribution from 
the pole of the Hurwitz zeta function aX 2z — 1 — m = \. All the poles are simple 
unless the pole of at 2 = 1 -|- y is a non-negative integer. Hence, the integral in 
fl98|l becomes 

f(") (-i;a,52) RiRes i+Res_,^m+ V (-1)^ V CH(2fc - 1 - m, a)6^^ (99) 
^=2 ^-1+2 ^ fc!r(-2) 

where the sum arises from the simple poles of 1/ sin(7r2;). Note that 

r(e) = ^-7 + 0(e) (100) 

and 

a(-A;,a) = -:^±i^, forfcGN, (101) 

where Bm{ci) are the Bernoulli polynomials and 7 is the Euler-Mascheroni constant. 
Hence, the residue at z = ^ gives 

Res ,=?^I^b. (102) 
^=2 m + 1 ^ ' 

Finally, for the residue at 2; = 1 + y, we note that 

CH{l + e,a) = ^-^{a) + 0{e), (103) 

where "^{a) = ^lnr(a) is the digamma function, see [25l p. 271]. The behaviour 
depends on whether m is even or odd. For odd m this is just another simple pole, 
and we obtain 



Y ' ' 4r(2+ 2. 



^2=1+ 

However, for even m there is a double pole, and we have to use 



^=^+^ z^i+fdzy 2; r(2+ l)r(-|)sin(7rz) ' ^ ' 

to obtain 



Res„.,™ = (-1) . 4„p^„)^fp,.M) (106) 

* ((^(Hi±l) _ ^(^) _ 2^(a) + 2 ln(6)) m(m + 2) - 4(1 + m)) . 
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We now use the same trick as in flHTI) to rewrite the regularized energy in fl92|) as 

En,)^o = -4F(^n-|;f,/) + (lor) 

The regularized energy is plotted in figure O As a consistency check it can be shown 
that Ef(^^)=Q{g = 0) = ^ as calculated in ( l90l) . It would be interesting to compare 
these results to other regularization methods. 



0.07- 
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0.5 0.75 
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Figure 2: The Dirac sea energy -E'/=o as a function of the coupling constant g. 



Finally, we address the problem of calculation the Dirac sea energy for fermions 
coupled to a i? = 1 background Skyrmion. In this case the Dirac Sea energy is given 
by 



(108) 



which is the sum of Egym and Easym for both parities. The energy of the "zero mode" 
Eq = ^ — g also needs to be taken into account, and we expect a similar picture as 
in [m. 
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Unfortunately, this is a much more comphcated situation, and zeta functions of 
this type have not been discussed in the hterature, to our knowledge. As a starting 
point, we could again perform a binomial expansion. We can then rewrite the energy 
Ef{^^)=^ as an infinite sum of zeta functions -F('^)(— a; IP'). Unfortunately, the last 
term in fllOSp leads to 6^ = —2g + g"^ which is negative for small and our formula 
no longer converges. It would be interesting to derive alternative expressions for 
these types of zeta function. 

5 Conclusion 

In this paper we consider the Dirac equation for fermions on chirally coupled 
to a spherically-symmetric background Skyrmion with topological charge one. The 
time-independent Dirac equation commutes with the grand spin and parity, and 
these symmetries allow us to reduce the Dirac equation to a system of four linear 
ODEs. Making use of the theory of Fuchsian differential equations, we derive the 
complete energy spectrum and the corresponding eigenfunctions which are given by 
polynomials. There is a positive parity state with energy 

and a negative parity state with energy 

EQ = l + g. 

Both states are generally non-degenerate. The energies 

E,%^(n) = \ ± ^{n + iy-2g + g^ 

and 

= 1 ± ^(n + lf + 2g + g^ 

all have degeneracy {n + 1)^, and correspond to states with parity (—1)*^ and 
(— l)*^"*"^, respectively. For G = 0, these energies were found in [22]. Finally, the 
energies 

E%^{n) = -\± ^{n + lf + g^ 

have degeneracy 2n{n + 2). The factor of 2 arises because the energy of these states 
is independent of parity. Furthermore, these states only occur for G > 0. 

For zero coupling {g = 0) the energy spectrum is E = ±{N + |) and the degen- 
eracy was found to be 

D = A{N + 1){N + 2), (109) 

in agreement with [27]. We also found explicit formulae for accidental degeneracies 
which occur for special values of the coupling constant g. 

The explicit formulae for the energy spectrum and its degeneracy enabled us 
to write down the zeta function related to the Dirac sea. For massive fermions on 
5*^, we were able to derive an asymptotic formula for a zeta function of generalized 
Epstein-Hurwitz type. The more interesting case of fermions coupled to Skyrmions 
on 5*^ leads to an interesting novel type of zeta function. However, we were unable 
to evaluate it using our current approach. This is an interesting topic for further 
study. 
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